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PREFACE 


This  report  was  prepared  as  part  of  a long-range  effort  to  tabu- 
late pertinent  climatological  and  topographic  information  for  the 
global  atmosphere  and  ocean  in  a model-compatible  format.  The  data 
are  used  in  the  work  of  the  Fiand/ARPA  Climate  Dynamics  Program,  spon- 
our*>d  by  the  Defense  Advanced  Research  Projects  Agency.  Related  tabu- 
lations are  those  of  ocean  depths  in  Rand  Report  R-1277-ARPA  and  of 
terrain  heights  in  R-127C-APPA.  Tabulations  of  global  atmospheric 
climatic  data  have  been  given  in  R-915-ARPA,  R-1029-ARPA,  R-1317-ARPA, 
and  R-1425-ARPA. 
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SUMMARY 


Climatological  monthly  ocean-surface  temperatures  obtained  from 
the  National  Center  for  Atmospheric  Research  and  from  Fleet  Numerical 
Weather  Central  are  merged  and  interpolated  onto  a 1°  global  grid. 
Monthly  distributions  of  the  main  ice  packs  of  the  Arctic  and  Antarctic 
are  digitized  from  Fleet  Weather  Facility  ice  charts  and  Navy  Atlases, 
and  are  incorporated  into  the  global  arrays.  Machine-analyzed  maps  of 
the  12  monthly  distributions  are  presented,  together  with  maps  and 
tabulations  of  averages  of  these  data  on  a global  grid  of  4°  latitude 
* 5°  longitude  for  the  months  of  February  and  August. 
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I . INTRODUCTION 


In  numerical  simulations  using  atmospheric  general  circulation 
models  (GCMs)  the  distributions  of  sea-surface  temperature  and  sea 
ice  are  required  as  input  boundary  conditions.  When  simulation  tiroes 
are  comparable  to  one  month  or  less,  it  is  usually  adequate  to  hold 
temperature  and  ice  distributions  fixed.  For  longer  time  scales,  as 
in  seasonal  and  interannual  climatic  simulations,  there  is  a need  for 
monthly  distributions  with  which  to  simulate  changing  boundary  condi- 
tions at  the  earth’s  surface.  With  the  development  of  oceanic  GCMs 
leading  toward  climatic  experiments  with  interactive  atmosphere/ ocean 
models,  there  is  the  additional  need  for  global  oceanic  data  with  which 
to  calibrate  and  test  both  the  ocean  models  and  the  coupled  models. 
Generally,  such  data  are  available  only  for  selected  ocean  basins  and 
are  not  in  a modol-compat ible  format. 

Recognizing  such  needs,  Washington  and  Thiel  (1970)  have  presented 
monthly  average  ocean-surface  temperatures  that  they  digitized  on  a 
2.5°  global  grid  from  a composite  of  sources.  Although  these  data  un- 
doubtedly constitute  the  most  complete  single  source  of  sea-surface 
temperatures  available,  they  do  not  extend  into  latitudes  higher  than 
approximately  60°N  and  S,  and  there  is  an  additional  need  for  ice-pack 
distributions.  The  purpose  of  the  work  described  here  is  to  revise  and 
extend  the  data  of  Washington  and  Thiel  to  the  entire  globe  on  a 1°  grid, 
and  to  include  in  the  tabulations  the  monthly  Jistr ibut ions  of  the  main 
ice  packs.  The  purpose  of  the  1°  grid  is  to  facilitate  application  to 
various  GCM  grids. 

Because  a tabulation  of  12  monthly  distributions  on  a 1°  global 
grid  is  too  large  for  convenient  publication,  wo  give  here  global  maps 
of  the  fine-mesh  data,  ogether  wi^h  tabulation  and  maps  of  averages 
of  these  data  on  a mesh  of  4°  latitude  “ 5"  longitude  for  the  months  of 
February  nr  August..  This  is  done  in  the  belief  that  these  maps  and 
tables  ma-  be  of  interest  to  the  meteorological  and  oceanographic  com- 
munities i.  general.  The  more  complete  data  on  the  1°  grid  are  avail- 
able on  magnetic  tape  upon  request. 


T.  SOURCES  AND  METHODS 


SEA  ICE 

For  the  Arctic  we  used  the  monthly  ice  charts  of  the  U.S.  Fleet 
Weather  Facility  (FWF).  These  charts,  which  give  distributions  of  ice 
greater  than  0,5  fractional  coverage,  are  in  general  agreement  with 
the  Hydrographic  Office  Atlas  (1S58)  but  are  less  detailed  in  that  they 
do  not  include  other  fractional  coverages.  The  advantage  of  the  FWF 
charts  is  that  they  apply  to  a larger  region  than  tha’  of  the  Hydro- 
graphic  Office  Atlas,  including  northern  portions  of  the  Pacific, 
Atlantic,  and  Hudson  Bay. 

The  rWF  charts  give  monthly  distributions  for  the  months  of  November 
through  March  and  semi-monthly  distributions  for  April  through  October. 
Rather  than  average  the  semi-monthly  distributions,  we  chose  the  ice 
charts  appropriate  to  the  lost  half  of  each  of  these  months  as  being 
sufficiently  representative.  This  introduces  a slight  bias  toward  more 
open  water  in  the  spring  and  toward  more  ice  in  the  fall,  but  has  the 
advantage  that  the  distributions  for  the  first  half  of  each  of  these 
months  can  be  included  later  without  requiring  modification  of  the 
existing  data. 

For  the  Antarctic,  we  used  the  monthly  mean  limits  of  the  ice  pack 
as  provided  in  the  Hydrographic  Office  (1957)  charts.  These  limits 
generally  correspond  to  ice  coverages  greater  than  0.5,  but  also  include 
some  c /ncentrations  within  the  boundary  of  the  main  pack  that  are  less 
than  0.5.  These  exceptions  occur  near  the  Antarctic  coast  during  the 
summer  months  of  December  through  February. 

The  ice  limits  were  digitized  by  hand  from  tht  r -,regoing  charts 
by  means  of  a piece-wise  linear  approximation  to  a continuous  line 
surrounding  the  pole.  Tne  boundaries  jl  disconnected  patches  of  ice 
were  also  included.  Wherever  the  ice  limit  Intersected  the  coast,  the 
line  was  continued  across  land  by  the  most  direct  route  to  the  nearest 
coastal  point  where  the  ice  limit  moved  out  to  sea  again.  From  input 
coordinates  of  the  ice  limits,  a computer  program  assigned  ice  to  those 
points  on  the  grid  that  were  completely  bounded  by  the  lines.  Particular 
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care  had  to  be  given  to  ihe  complicated  summer  distributions  in  the 
Canadian  Arctic.  We  believe  that  the  results  are  accurate  and  realistic 
when  a realistic  land  distribution  is  used  to  suppress  lrnd  points. 

SF- SURFACE  TEMPERATURES 

Global  monthly  sea-surface  temperature  data  were  obtained  on  mag- 
netic tape  from  the  National  Center  for  Atmospheric  Research  (NCAR) , 
and  are  as  documented  by  Washington  and  Thiel  (1970).  Valid  data  points 
were  represented  by  Centigrade  temperatures  on  the  NCAR  2.5°  global  grid. 
Missing  d^ca  points  in  latitudes  higher  than  approximately  60°N  and  S, 
as  well  as  all  land  points,  were  designated  by  99°. 

Our  first  task  was  to  interpolate  or  extrapolate  the  NCAR  data 
into  land  areas  and  high  latitudes  in  order  to  facilitate  further  inter- 
polation onto  a 1°  grid.  The  Hydrographic  Office  Atlas  (±957)  reports 
that  water  temperatures  beneath  the  Antarctic  ice  pack  are  typically 
29°F  (-1.7°C)  throughout  the  year.  Havin'  digitized  the  monthly  ice 
limits,  we  assigned  a uniform  temperature  of  -1.7°C  to  all  NCAR  mesh 
points  lying  within  areas  covered  by  sea  ice,  according  to  month.  The 
remaining  points  with  missing  data  were  filled  in  by  a form  of  bilinear 
interpolation,  defined  as  the  average  of  east-west  and  north-south  linear 
interpolations  between  the  nearest  valid  data  points.  The  linear  inter- 
polations were  weighted  according  to  distance  from  the  valid  datum.  In 
some  instances,  only  a single  interpolation  in  one  direction  was  possible. 
Two-dimensional  Bessel  interpolation  was  then  used  (Slnffensen,  1950, 
p.  210)  to  fill  in  a 1°  global  mesh,  giving  temperatures  at  all  inter- 
sections of  integer  latitude  and  longitude.  At  this  point,  however,  we 
did  not  regard  the  results  north  of  approximately  bO°N  as  being  valid. 

A better  representation  was  needed  for  the  Arctic  and  adjacent  seas. 

The  second  major  source  of  monthly  sea-surface  temperatures  was  the 
northern- hemisphere  climatological  data  that  we  obtained,  also  on  mag- 
netic tape,  from  U.S.  Fleet  Numerical  VJeather  Centra'  (F.'WC)  . These 
data,  in  current  use  by  FNVC  as  of  late  1972,  were  defined  on  the  Navy 
grid,  which  is  a 125  * 125  rectangular  grid  superimposed  on  a poiar 
sterographic  projection,  true  at  60°N.  The  equator  is  an  inscribed 
circle,  meaning  that  the  rectangular  mesh  extends  somewhat  south  of 
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the  equator  in  the  vicinity  of  four  meridians.  These  meridians  are 
at  125CW,  35°W,  55°E,  and  145°E,  where  grid  points  extend  to  approxi- 
mately lb°S.  7c lues  had  already  been  interpolated  onto  land  points 
by  FNWC,  so  no  further  interpolation  of  this  type  was  required.  Again, 
Bessel  interpolation  was  usee'  to  transom  the  data  onto  the  1°  grid. 

A decision  was  required  at  this  stage  concerning  where  to  merge 
the  two  sets  of  data:  (a)  ir  high  northern  latitudes  where  the  NCAR 

data  are  missing  or  \,o)  near  the  equator  and  low  «■  mthern  latitudes 
where  the  FNWC  data  terminate.  There  are  other  possible  approaches, 
but  we  chose  (b)  because  the  FNWC  data  appeared  to  give  better  resolu- 
tion of  the  cold  equatorial  region  in  the  eastern  Pacific  in  January 
and  February.  The  method  of  merging  employed  weighted  averages  of  the 
NCAR  and  FNWC  data  within  a 10°  latitude  belt  bounded  to  the  south  by 
the  limits  of  the  Navy  grid.  The  weighting  was  linear  with  latitude. 

As  a final  step,  we  applied  a nine-point  smoothing  operator  (Shapiro, 
1970)  to  the  resulting  12  global  arrays  to  remove  a small  amount  of 
noise  introduced  by  the  interpolation  and  merging  procedures. 
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III.  RESULTS 


By  the  methods  described  in  Sec.  II  we  obtained  12  monthly  18i 
x 360  arrays  with  temperatures  stored  at  every  integer  value  of  lati- 
tude and  longitude.  Temperatures  for  water  that  is  covered  by  ice 
are  given  in  degrees  Kelvin,  for  which  we  have  used  the  transformation 
273. 1°K  = 0.0°C.  The  remainder  are  in  degrees  Centigrade.  This  was 
merely  a cor. lenient  way  to  distinguish  ice-covered  wate’*  from  open 
water  in  our  data.  As  described  in  Sec.  II,  the  ice-pack  distributions 
cross  land  in  the  basic  data  files,  but  the  results  are  ft,  . id  to  be 
realistic  when  land  points  are  masked  out  by  a realistic  land  distribu- 
tion. Machine-analyzed  maps  of  the  1°  mesh  data  are  given  in  Figs.  1 
through  12,  in  which  ice-pack  distributions  are  designated  by  broken 
hash  marks.  The  land  outlines  are  based  on  the  data  of  Smith  et  al. 
(1966),  as  corrected  and  interpolated  onto  a 1°  global  mesh  by  Gates 
end  Nelson  (1973).  Gall's  stereographic  projection  (Bartholomew,  1966), 
which  is  true  at  45°N  and  S,  was  used  for  the  figures. 

Concerning  resolution,  we  emphasize  that  use  of  a 1°  global  mesh 
is  for  convenience  in  interpolating  onto  various  other  grids.  No 
accuracy  greater  than  that  implied  by  the  NCAR  2.5°  global  mesh,  or 
the  FNWC  125  x 125  mesh,  can  be  claimed  for  the  temperature  data  pre- 
sented here.  However,  we  believe  that  the  accuracy  of  the  main  ice-pack 
distributions  is  consistent  with  that  implied  by  a 1°  mesh. 

The  1°  mesh  data  have  been  area-averaged  onto  the  46  * 72  global 
grid  (4°  latitude,  5°  longitude)  that  is  currently  used  in  general 
circulation  modeling  at  The  Rand  Corporation,  UCLA,  and  the  Goddard 
Institute  for  Space  Studies.  The  results  for  the  months  of  February 
and  August  are  given  in  Figs.  13  and  14,  and  in  Tables  1 and  2.  The 
land  distributions  of  Figs.  13  and  14,  which  correspond  to  the  asterisks 
of  Tables  1 and  2,  were  objectively  determined  from  tne  1°  data  of  Gates 
and  Nelson  (1973).  This  is  the  land  outline  in  current  use  with  the 
Rand  atmospheric  GCM. 
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IV.  RECOMMENDED  FUTURE  OCEANIC  DATA  ACQUISITI ON<- 


As  oceanic  GCMs  become  more  sophisticated,  e.g.,  have  increased 
resolution  or  include  ice-pack  models,  there  will  be  increasing  re- 
quirements for  more  detailed  oceanic  data  in  a model-compatible  format. 
Some  suggested  categories  for  needed  improvements  are  described  below. 

SEA  ICE 

A high-priority  item  in  this  category  is  to  include  the  remaining 
semi-monthly  distributions  for  the  months  of  April  through  October  in 
the  Arctic.  This  would  give  an  improved  representation  of  the  more 
rapidly  changing  ice  distributions  that  occur  in  the  warmer  months. 
Similar  details  should  be  added  tor  the  Antarctic  when  available.  In 
addition,  more  detailed  fractional  coverages  of  ice  (e.g.,  <0.1,  0.1 
to  0.5,  and  so  forth)  will  probably  be  required.  The  degree  of  cover- 
age and  the  presence  or  absence  of  openings  and  leads  are  important  in 
computing  atmosphere/ocean  heat  transfer. 

SEA-SURFACE  TEMPERATURES 

Anticipating  the  nepd  for  increased  horizontal  resolution  of  global 
ocean  models,  we  will  want  to  add  n.ore  detailed  temperature  distribu- 
tions consistent  with  l°-mesh  accuracy.  An  example  is  the  monthly  dis- 
tributions of  sea-surface  temperature  for  the  eastern  Pacific  compiled 
by  Wyrtki  (1964) . 

Measurements  of  sea-surface  temperature  by  earth  satellite  repre- 
sent a promising  way  to  obtain  globax  synoptic  data.  As  methods  are 
developed  for  increasing  the  accuracy  of  these  measurements  and  filter- 
ing out  the  effects  of  cloud  cover,  the  addition  of  sich  data  to  a 
digitized  climatological  data  base  should  be  valuable,  particularly 
in  sparse-data  areas. 

TEMPERATURE  ANr  SALINITY  WITH  DEPTH 

A compilation  of  global  sea-surface  temperatures  iu  only  a begin- 
ning task  in  the  process  of  building  an  adequate  climatological  data  base 
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for  the  global  ocean.  The  density  and  freezing  point  of  sea  water 
depend  on  salinity  as  well  as  temperature.  Perhaps  the  highest  priority 
in  future  work  should  be  given  to  obtaining  at  .least  seasonal  averages 
of  the  global  distributions  of  both  temperature  and  salinity  at  standard 
depths  for  tne  upper  ocean,  say,  the  upper  200  to  300  m.  Much  of  the 
required  raw  data  are  available  from  the  National  Oceanographic  Data 
Center  (NODC,  1969),  but  these  will  have  to  be  analyzed,  and  there  are 
appreciable  areas  in  which  data  are  sparse. 


ig.  1 January  average  ocean-surface  temperatures  (°C)  and  ice-pack  distributions 

(broken  hash  marks).  Solid  isolines  denote  0°,  2°,  1°,  6°C  ...  . Dotted  isolines 
in  low  latitudes  denote  25°,  27°,  and  2JTC. 
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Fig.  5 May  average  ocean-surface  temperatures  (°C)  and  ice-pack  distributions 

(broken  hash  marks).  Solid  isolines  denote  0°,  2°,  4°,  6°C  ....  Dotted  isolines 
in  low  latitudes  denote  25°,  27°,  and  29°C. 
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Fig.  10 October  average  ocean-surface  ♦'"nperatures  (°C)  and  ice-pack  distributions 

(broken  hash  marks).  Solid  isolines  denote  0°,  2°,  4°,  6°C  ...  . Dotted  isolines 
in  low  latitudes  denote  25°,  27°,  and  2b°C. 
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11 November  average  ocean-surface  temperatures  (°C)  and  ice-pack  distributions 

(broken  hash  marks).  Solid  isolines  denote  0°,  2°,  4°,  6°C  ...  . Dotted  isolines 
in  low  latitudes  denote  25' , 27°,  and  29°C. 
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roken  hash  marks).  Solid  isolines  denote  0°,  2°,  4°.  6°C  . . 

in  low  latitudes  denote  25°,  27°,  and  29°C. 


F;g.  13 — February  average  ocean-surface  temperatures  (°C)  and  ice-pack  distributions 
(broken  hash  marks),  averaged  from  basic  data  for  mesh  of  4°  lat  X S3  long.  Solid 
isolines  denote  0°,  2°,  4°C  ...  . Dotted  isolines  denote  25°,  27°,  and  29°C. 


ig.  14 August  average  ocean-surface  temperatures  (°C)  and  ice-pack  distributions 

(broken  hash  marks),  averaged  from  basic  data  for  mesh  of  4°  lat  X 5°  long.  Solid 
isolines  denote  0°,  2°,  4°C  ....  Dotted  isolines  denote  25°,  27°,  and  29°C. 
Grid-point  data  are  given  in  Table  2. 
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Table  1 — Continued 
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AUGUST  AVERAGE  UCEAN-SUR7ACE  TEMPERATURE^ . OPEN-WATER  TEMPERATURES  IN  °i 
WATER  TEMPERATURES  BENEATH  SEA  ICE  IN  °K.  ASTERISKS  DENOTE  LAND  POINTS 
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